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Chiral effective field theory (xEFT) complements numerical simulations of quantum chromody- 
namics (QCD) on a space-time lattice. It provides a model-independent formalism for connecting 
lattice simulation results at finite volume and a variety of quark masses to the physical world. The 
asymptotic nature of the chiral expansion places the focus on the first few terms of the expansion. 
Thus, knowledge of the power-counting regime (PCR) of xEFT, where higher-order terms of the ex- 
pansion may be regarded as negligible, is as important as knowledge of the expansion itself. Through 
the consideration of a variety of renormalization schemes and associated parameters, techniques to 
identify the PCR where results are independent of the renormalization scheme are established. The 
nucleon mass is considered as a benchmark for illustrating this general approach. Because the PCR 
is small, the numerical simulation results are also examined to search for the possible presence of 
an intrinsic scale which may be used in a nonperturbative manner to describe lattice simulation 
results outside of the PCR. Positive results that improve on the current optimistic application of 
chiral perturbation theory (\PT) beyond the PCR are reported. 

PACS numbers: 12.39.Fe 12.38.Aw 12.38.Gc 14.20.Dh 



I. INTRODUCTION 

The low energy chiral effective field theory (xEFT) 
of quantum chromodynamics (QCD) provides a model- 
independent approach for understanding the conse- 
quences of dynamical chiral-symmetry breaking in the 
chiral properties of hadrons. Non-analytic contributions 
in the quark mass are generated by the pseudo-Goldstone 
meson dressings of hadrons through meson-loop integrals. 
Chiral perturbation theory (%PT) provides a formal ap- 
proach to counting the powers of low energy momenta 
and quark masses such that an ordered expansion in pow- 
ers of the quark mass m q oc m 2 is constructed. x?T indi- 
cates that, in general, the most singular nonanalytic con- 
tributions to hadron properties lie in the one-loop 'meson 
cloud' of the hadron. For example, the leading nonana- 
lytic behavior of a baryon mass is proportional to rrCq 7, 
or rrv^. More generally, baryon masses can be written as 
an ordered expansion in quark mass or m 2 . 

To establish a model-independent framework in xPT, 
the expansion must display the properties of a convergent 
series for the terms considered. There is a power-counting 
regime (PCR) where the quark mass is small, and higher- 
order terms in the expansion are negligible beyond the 
order calculated. Within the PCR, the truncation of the 
chiral expansion is reliable to a prescribed precision. 

The asymptotic nature of the chiral expansion places 
the focus on the first few terms of the expansion. A 
survey of the literature for the baryon sector of xEFT 
illustrates the rarity of calculations beyond one-loop 1 
3], and there are no two- loop calculations which incor- 
porate the effects of placing a baryon in a finite volume. 
With only a few terms of the expansion known for cer- 
tain, knowledge of the PCR of xEFT is as important as 
knowledge of the expansion itself. It is within the PCR 
that higher-order terms of the expansion may be regarded 



as negligible. 

Numerical simulations of QCD on a space-time lattice 
are complemented by %EFT through the provision of a 
model-independent formalism for connecting lattice sim- 
ulation results to the physical world. Simulations at finite 
volume and a variety of quark masses are related to the 
infinite volume and physical quark masses through this 
formalism. However, the formalism is accurate only if 
one works within the PCR of the truncated expansion. 
Present practice in the field is best described as opti- 
mistic. Truncated expansions are regularly applied to a 
wide range of quark (or pion) masses with little regard 
to a rigorous determination of the PCR. 

When considering nucleons, there is some evidence 
that the PCR may be quite small; constrained by 
to tt S 200 MeV at 1% accuracy at the chiral order 
0{m\ log m^) H 01 ■ Th is estimate of the PCR of xPT 
was identified using specific finite-range regularization 
(FRR) techniques to analyze lattice QCD data. Using 
FRR, the regime is manifest when the quark-mass depen- 
dence of the nucleon mass is independent of the renor- 
malization scheme parameter. 

A chief focus of this paper is to establish a rigorous 
approach to determining the PCR of a truncated chiral 
expansion quantitatively. Through the consideration of 
a variety of renormalization schemes and associated pa- 
rameters, new techniques to identify the PCR are estab- 
lished. The PCR is the regime where results are scheme 
independent. The nucleon mass is considered as a bench- 
mark for illustrating this general approach. Here, the 
chiral expansion is examined, focusing on the individual 
low energy coefficients of the chiral expansion. This ap- 
proach provides a determination of the PCR for a trun- 
cated expansion in xEFT. As discussed in detail in the 
following section, the PCR is indeed small for the nucleon 
mass. Other observables are expected to show a similar if 
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not smaller PCR. Thus most of today's lattice simulation 
results lie outside the PCR, and the truncated chiral ex- 
pansions have been used to extrapolate from outside the 
PCR. The low energy coefficients determined by apply- 
ing the truncated expansion outside the PCR will take 
on unphysical values, as they accommodate important 
but otherwise missing contributions from non-negligible 
higher-order terms. 

While continued advances in numerical simulations of 
lattice QCD will be vital to some extent in resolving this 
problem, the physical value of the strange-quark mass 
presents a challenge that will not diminish with super- 
computing advances. If one were to include the effects of 
kaons, vital to understanding strangeness in the nucleon 
for example, then one must either calculate to signifi- 
cantly higher order in the expansion of x^T, or develop 
new non-perturbative approaches which utilize the non- 
perturbative information expressed in the lattice simula- 
tion results. Since the former is likely to be compromised 
by the asymptotic nature of the expansion, attention is 
given to the latter approach. 

Thus the second focus of this paper is to examine 
the numerical simulation results, to identify the possi- 
ble presence of an intrinsic scale. This may then be 
used to address lattice simulation results outside of the 
PCR in a non-perturbative manner. Of course, the non- 
perturbative formalism must incorporate the exact per- 
turbative results of xPT in the PCR. Positive results are 
reported that improve on the current optimistic applica- 
tion of xPT outside of the PCR. 

The outline of the presentation is as follows. Sec. |TT] re- 
views chiral effective field theory and the process of reg- 
ularization and renormalization. The adoption of finite- 
range regularization (FRR) provides a wide range of 
schemes and scales, which overlap with the more pop- 
ular massless renormalization schemes as the finite-range 
regulator parameter is taken to infinity. Sec. IIIII inves- 
tigates FRR in the context of a particular model. By 
generating a set of pseudodata and analysing it with a 
variety of renormalization schemes, a robust method for 
determining the PCR is obtained, along with an optimal 
renormalization scale to use beyond the PCR. Finally, 
Sec. IIVI includes the analysis of three sets of lattice re- 
sults for the nucleon mass, utilizing the tools developed 
in the previous section. Conclusions are summarized in 
Sec. El 



II. EFFECTIVE FIELD THEORY 

This section begins by briefly reviewing the process of 
regularization and renormalization in finite-range regu- 
larized (FRR) chiral effective field theory, providing a 
range of renormalization schemes and scales. A central 
focus is to search for the dependence of physical results 
on the scheme and associated scales, as these will be an 
indication that one is applying the chiral expansion out- 
side the PCR. 



The focus is to establish techniques that provide a 
quantitative test of whether a given range of m n lies 
within the PCR. This is achieved through an examina- 
tion of the flow of the low energy coefficients as a function 
of the renormalization scheme parameter(s). A negligi- 
ble dependence would confirm that the pion-mass range 
is within the PCR. On the other hand, the properties of 
the flow will be used to identify a preferred regularization 
scheme in a non-perturbative sense that best describes 
the results beyond the PCR. 

A. Renormalization in FRR xEFT 

Using the standard Gell-Mann-Oakes-Renner relation 
connecting quark and pion masses, m q oc toJ [6], the 
formal chiral expansion of the nucleon can be written 
as a polynomial expansion in plus the meson-loop 
integral contributions: 

M N = {a + a 2 ml + a 4 m* + O(m^)} 

+ E w + E A + E tad . (1) 

The pion cloud corrections are considered in the heavy- 
baryon limit, with loop integrals, E^r, Ea and E ta( ;, cor- 
responding to Figures [1] through 02 respectively. The 
coefficients of the analytic polynomial, contained in 
brackets { } in Eq. (fTJ), are related to the low energy 
constants of xPT. In this investigation, they will be de- 
termined by fitting to lattice QCD data. These coef- 
ficients will be referred to as the residual series coeffi- 
cients. These bare coefficients undergo renormalization 
due to contributions from the loop integrals E^v, Ea and 

Efad- 

Under the most general considerations, each loop in- 
tegral, when evaluated, produces an analytic polynomial 
in and nonanalytic terms: 

E W = b$ + b^ml + XN ml + b?mi + 0(rnl) , (2) 
E A = b£ + b£ml + b£mi 

+ ^ K XAmilog^ + 0(ml), (3) 

E tad = &|V4 + &i'rr4+xXl°g— + 0(m«).(4) 

A* 

Here A is the delta-nucleon mass splitting in the chiral 
limit, taken to be 292 MeV. xn, Xa and x't denote the 
model independent chiral coefficients of the terms that 
are nonanalytic in the quark mass. The 6.; coefficients 
are renormalization scheme dependent as are the co- 
efficients. It can be noted that the tadpole loop contri- 
bution T<tad does not produce a 6q term because it enters 
with a leading factor of mj, as discussed in Section Hi Bl 
The primes on the coefficients b\ and x't ncrc simply in- 
dicate that they will be used later in a slightly different 
form. 

The process of renormalization in FRR %EFT proceeds 
by combining the renormalization-scheme dependent co- 
efficients to provide the physical low energy coefficients, 
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FIG. 1: The pion loop contribution to the self-energy of the 
nucleon, providing the leading nonanalytic contribution to the 
nucleon mass. All charge conserving transitions are implicit. 
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FIG. 2: The pion loop contribution to the self-energy of 
the nucleon allowing transitions to the near-by and strongly- 
coupled decuplet baryons. 
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FIG. 3: Tadpole contributions to the nucleon self energy. 

which are denoted as Cj. Thus, the nucleon mass expan- 
sion takes on the standard form: 

M N = c + c 2 ml + XNml + c±m% 

+ (-^XA + caxtjmJlog^+OK). (5) 

By comparing Eqs. (JIJ through ([5]), the following renor- 
malization procedure is obtained: 

co - a + b$ + b£, (6) 
c 2 = a 2 + b$ + b$ + bt, (7) 
c 4 = a 4 + bf + bf + b{ , etc. (8) 

The coefficients Ci are scheme independent quantities, 
and this property will be demonstrated when determined 
within the PCR. The value of Co is the nucleon mass in 
the chiral limit (m^ = 0), and c 2 is related to the so- 
called sigma term of explicit chiral symmetry breaking 
The nonanalytic terms and m\ logm w //x have 
known, constant coefficients denoted by xn, Xa and Xt- 
The value of c 4 is scale dependent, such that the total m% 
term in Eq. ([5]), including the logarithm, is independent 
of the scale /i. It can be noted that the nucleon mass 
itself is completely independent of the choice of \i. For 
the numerical analysis, /i is set equal to 1 GeV. 

Of course, EFT loop calculations are commonly di- 
vergent without some regularization method. Since the 
effective field theory is only applicable for low energies, 
hard momenta contributions to loop calculations may be 
eliminated. However, the traditional schemes including 
dimensional regularization (DR) often do not involve an 



explicit scale dependence when evaluating loop diagrams. 
Without any momentum cutoff, the bi coefficients from 
each loop integral become either infinite or vanish, and 
the Ci coefficients from Eq. ([5]) undergo an infinite renor- 
malization or none at all: 

co = a - oo , (9) 
c 2 = a 2 + oo, (10) 
c 4 = a 4 + , etc. (11) 

Since the c 2 ; coefficients are finite after renormalization, 
the at coefficients must have been infinite, with the oppo- 
site sign of the bi coefficients. As emphasized above, both 
the ai and bi coefficients are scheme dependent. The in- 
finities are absorbed in constructing the Cj coefficients 
and thus subtracted from the chiral expansion. This 
minimal subtraction scheme with no explicit scale de- 
pendence makes DR quite suitable for elementary fields, 
where the absence of new degrees of freedom at higher 
energies is assumed. However, for EFTs there is an en- 
ergy scale beyond which the effective fields are no longer 
the relevant degrees of freedom. When one integrates 
loop contributions over this high energy domain, there is 
no guarantee that one can efficiently subtract the model- 
dependent, ultraviolet physics with a finite number of 
counter-terms (unless in the PCR). As a result, the chi- 
ral expansion typically only shows reliable convergence 
properties over a narrow range of pion mass. 

Indeed this problem of beginning with rapidly varying 
loop contributions, which must then be removed with a 
finite number of counter-terms, can easily be overcome. 
The hard momentum contributions to the meson-loop di- 
agrams can be suppressed via the introduction of a regu- 
lator. As such, the coefficients of the residual expansion 
are likely to be smaller, and the utility of the expansion 
has the potential to apply to a broader range of quark 
or pion masses. The introduction of a regulator acts to 
resum the chiral expansion, with loop integrals having 
the general properties described in Eqs. © through ((U). 

The resummation of the chiral series through the intro- 
duction of a regulator (or similar variant) has been stud- 
ied in various instances The method consists 
of inserting a regulator function u(k 2 ) into the integrand 
of the meson-loop integrals. The regulator can take any 
form, so long as it is normalized to 1, and approaches 
sufficiently fast to ensure convergence of the loop. Unlike 
DR, this method involves an explicit momentum cutoff 
scale, A. The chiral expansion can now be written in 
terms of this cutoff scale: 

Mn = {a^ + 4ml+a^ml + 0{ml)} 

+ V N {ml , A) + E A ( TO 2 > A ) + i: tad (ml , A). (12) 

The superscript A denotes the scheme dependence of the 
coefficients. The loop integrals are functions of the 
scale A and also m\. 

Through the introduction of the regulator, the loop in- 
tegrals are now low energy contributions, significant for 
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small m 2 and becoming negligible as m 2 . becomes large. 
The scheme dependent af' coefficients undergo a renor- 
malization, as before, via their combination with the bf 
coefficients, whose scheme dependence is now explicit, 
reflecting the regularization of the loop integrals: 
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(13) 
(14) 
(15) 



Dimensional analysis reveals that the coefficients 6; are 
proportional to A^ 3 ~ l \ Thus it can be realized that as 
the cutoff scale A goes to infinity the FRR expansion 
reduces to that of Eq. © via Eqs. © through (HI]). At 
any finite A, a partial resummation of higher-order terms 
is introduced. 

Previous studies indicate that extrapolation results 
show very little sensitivity to the precise functional form 
of the regulator [lij]. In this investigation, the family 
of smoothly attenuating dipole regulators will be con- 
sidered. The general n— tuple dipole function takes the 
following form, for a cutoff scale of A: 



k 2n 

^)=( 1 + A^ 



(16) 



The standard dipole is recovered for n = 1. The cases 
n = 2, 3 are the 'double-' and 'triple-dipole' regulators, 
respectively. In the following, u{k 2 ) is used to denote 
one of these regulators. This functional form allows one 
to interpolate between the dipole regulator and the step 
function (which corresponds ton-> oo). 

In a study by Bernard et al. [l5[, it was suggested 
that only a sharp cutoff FRR scheme is consistent with 
chiral symmetry. Djukanovic et al. [l6j have demon- 
strated more general functional forms can be generated 
by proposing a scheme in which the regulator function is 
interpreted as a modification to the propagators of the 
theory, obtained from a new chiral symmetry-preserving 
Lagrangian. Higher-derivative coupling terms are built 
into the Lagrangian in order to produce a regulator from 
the Feynman Rules in a symmetry preserving manner. 

The regulators used in the present investigation are 
introduced in a less systematic fashion, such that chiral 
symmetry is not automatically preserved to the order cal- 
culated. The higher derivative couplings of the regulator 
induces scheme-dependent nonanalytic terms. To main- 
tain chiral symmetry, one must introduce the necessary 
vertex corrections. 

Alternatively, one can choose the regulator judiciously 
such that any extra scheme-dependent nonanalytic are 
removed to any chosen order. For example, the n— tuple 
dipole regulators generate extra nonanalytic terms in the 
chiral expansion of Eq. (|S|) at higher chiral orders. For a 
dipole regulator, regulator-dependent nonanalytic terms 



occur at odd powers of m^, beginning at O(m^) [30]. 
In the case of the double dipole, the nonanalytic terms 
begin at 0(m£), and for the triple dipole the nonanalytic 
terms begin only at O(m^). 

In a final observation, it is essential to note that the 
degrees of freedom present in the residual series coeffi- 
cients, a^, are sufficient to eliminate any dependence on 
the regulator parameter, A, to the order of the chiral ex- 
pansion calculated: in this case 0(m^.). By definition, 
higher order terms in the FRR expansion are negligible 
in the PCR, and therefore FRR xEFT is mathematically 
equivalent to xPT in the PCR. Any differences observed 
in results obtained at the same chiral order but with dif- 
ferent regularization schemes are a direct result of con- 
sidering data that lie outside the PCR (provided that the 
regulator A is not chosen too small such as to introduce 
an unphysical low energy scale). 



B. Loop Integrals and Definitions 

The leading order loop integral contributions to the 
nucleon mass, corresponding to the diagrams in Figures 
[T] through [3J can be simplified to a convenient form by 
taking the heavy-baryon limit, and performing the pole 
integration for fco- Renormalization, as outlined above, is 
achieved by subtracting the relevant terms in the Taylor 
expansion of the loop integrals and absorbing them into 
the corresponding low energy constants, Cf. 
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(18) 
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"'- -^XA"4log— (20) 



J tad 



1 



3 2u 2 {k 2 ) A;t 



c 2 ml ( 64 'VJ + Xtml log • 



C2«4<5- tad . 



(21) 

- + 0(ml)) (22) 
(23) 



These integrals are expressed in terms of the pion energy, 
u>(k) = yjk 2 + m 2 . The tilde (~) denotes that the inte- 
grals are written out in renormalized form to chiral order 



[30] While scheme-dependent, it is significant to note that with a 
dipole regulator, A = 0.8 GeV, the coefficient of the induced 
miterm compares favorably with the two-loop calculation 0- 

lira. 
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O(m^). As the bi coefficients are regulator and scale de- 
pendent, this subtraction removes this dependence. The 
coefficients ao and ai of the analytic terms in the chiral 
expansion in Eq. ([5]) are now automatically the renor- 
malized coefficients cq and C2- This is because the bo 
and 62 terms in Eqs. (fT3|) and (fT4|) are removed in the 
subtraction. Note also that the tadpole amplitude in 
Eqs. ([25]) and (|23p contains the renormalized C2 in its 
coefficient. The interaction vertex in this diagram arises 
from expanding out the pion field in the leading quark 
mass insertion. 

The constant coefficients xn, Xa and \t for each in- 
tegral are defined in terms of the pion decay constant, 
which is taken to be / w = 93 MeV, and the axial cou- 
pling parameters D, F and C which couple the baryons 
to the pion field. The phenomenological values for these 
couplings are used, applying the SU(6) flavour-symmetry 

relations [H to y ield C = ~ 2£) ' F 
value D = 0.76: 



volume V = L x L y L z [22 



|£) and the 
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XA 

Xt 



c 2 , 



327T/2 9 

3 
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(24) 
(25) 
(26) 



With the renormalized integrals specified, the FRR 
modified version of the chiral expansion in Eq. ([5]) takes 
the form: 



M 



N 



c + c 2 m 2 w (l + a tad ) + a%m% + Y, N + t A . (27) 



The a% ter m is left in unrenormalized form for simplicity. 
Indeed, the 64 can be evaluated by expanding out corre- 
sponding loop integrals, such as in Ref. [l2j. However, 
the focus here is on the behaviour of cq and c%. 

Since the results of lattice simulations reflect the pres- 
ence of discrete momentum values associated with the fi- 
nite volume of the lattices, the formalism must also take 
into account these finite volume effects. In order to ac- 
commodate the effect of the finite volume, the continuous 
loop integrals occurring in the meson loop calculations in 
infinite volume are transformed into a sum over discrete 
momentum values. The difference between a loop sum 
and its corresponding loop integral is the finite volume 
correction, which should vanish for all integrals as m T L 
becomes large [2ll ] . 

While Eq. (|27]l is useful in describing the pion mass 
evolution of the nucleon mass, for the consideration of 
lattice QCD results, one also needs to incorporate cor- 
rections to allow for the finite- volume nature of the nu- 
merical simulations. As the pion is the lightest degree of 
freedom in the system, it is the leading order pion loop 
effects that are most sensitive to the periodic bound- 
ary conditions. The corrections can be determined by 
considering the transformation of each loop integral in 
Eqs. ([T7|) . (HHJ) and (|2"Tj) . into a discrete sum for lattice 



d^fc 



(2tt)- 



L x L y L z k . ^ 



E 



(28) 



Each momentum component is quantized in units of 
2tt/L, that is fej = rii2ir/L for integers rij. The finite 
volume correction S F can be written as the difference 
between the finite sum and the integral: 



S FVC = 



Xi 
2tt 2 



J2li(k,ml,A)- d 3 k Ii(k,ml,A) 



(29) 

where i = N,A, and the integrands are denoted 
h(k, to 2 , A). By adding the relevant finite volume correc- 
tion (FVC) to each loop contribution, the finite volume 
nucleon mass can be parameterized: 



M V N 



c + 



+ S N VC ) 



&tad) 
V (E A 



a 4 m„ 

cFVC' 



'A 



(30) 



It is also anticipated that the FVC are independent of 
the regularization scale A in this domain. In Figures [4] 
and [5j the scale dependence of the finite- volume correc- 
tions is shown for a dipole regulator and a 2.9 fm box 
(the same box size used for the PACS-CS data [Hj]). It 
is notable that choosing A too small suppresses the very 
infrared physics that one is trying to describe, and there- 
fore it is sensible to be cautious by not selecting a A 
that is too low. Figures [B] and [7] show the behaviour of 
the FVC for a 4.0 fm box, and the corrections are much 
smaller as expected. 

For large A the results saturate to a fixed result. For 
the light pion masses, provided A > 0.8 GeV, the esti- 
mated finite volume corrections are stable. The asymp- 
totic result is used, which has been demonstrated to be 
successful in previous studies [13]. Numerically, this is 
achieved by evaluating the finite volume corrections with 
a parameter, A' = 2.0 GeV, 6f vc = 5f vc {A'). It should 
be noted that this is equivalent to the more algebraic 
approach outlined in Ref. (2lj . 



III. INTRINSIC SCALE: PSEUDODATA 

This xEFT extrapolation scheme to order 
0(m\ logm-n-) will be used in conjunction with lat- 
tice QCD data from JLQCD, PACS-CS and CP-PACS 
to predict the nucleon mass for any value of mj. 
The lattice data used in this analysis can be used to 
extrapolate Afjv to the physical point by taking into 
account the relevant curvature from the loop integrals 
in Eqs. (fT5|) , (|20p and (f2"2"j) . As an example, a regulator 
value of A = 1.0 GeV was chosen for Figures [8] through 
[TOl where the finite volume corrected EFT appears 
concordant with previous QCDSF-UKQCD results [24j. 
If the regulator is changed away from the choice A = 1.0 
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FIG. 4: (color online). 
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FIG. 6: (color online). 
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FIG. 7: (color online). 




original data 

A = 1.0 GeV, dipole, fin.vol. 
A = 1.0 GeV, dipole, inf.vol. 
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FIG. 8: (color online). Example dipole extrapolation based on 
JLQCD data box size: 1.9 fm. 
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FIG. 9: (color online). Example dipole extrapolation based on 
PACS-CS data [H, box size: 2.9 fm. 



GeV, the extrapolation curve also changes. This signifies 
a scheme dependence in the result due to using lattice 
QCD data beyond the PCR. 

To demonstrate this, consider the infinite volume ex- 
trapolation of the CP-PACS data. The extrapolation is 
achieved by subtracting the finite volume loop integral 
contributions defined in Eqs. (|T7|) . (|T^|) and (JUJ) from 
each data point and then fitting the result to obtain the 
coefficients cq, and as shown in Eq. PT|). The in- 
finite volume loop integrals are then added back at any 
desired value of m\. 

Figure [TT] shows that the curves overlap exactly when 
m\ is large, where the lattice data reside. They diverge 
as the chiral regime is approached. This section addresses 
this problem in detail. 

A particular regularization scale is selected and a dense 
and precise data set is generated, which smoothly con- 
nects with state of the art lattice simulation results. If all 
the data considered lie within the PCR then the choice 
of regulator parameter is irrelevant, and the FRR chiral 
expansion is mathematically equivalent to scale-invariant 
renormalzation schemes including DR. However, the pur- 
pose here is to consider an insightful scenario, whereby 
a set of ideal pseudodata with known low energy coeffi- 
cients is produced. This scenario will form the basis of 
the investigation of the PCR, and ultimately the possi- 
ble existence of an intrinsic scale hidden within the actual 
lattice QCD data. 

The pseudodata are produced by performing a finite 
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FIG. 10: (color online). Example dipole extrapolation based on 
CP-PACS data box size: 2.3 - 2.8 fm. 
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FIG. 11: (color online). Close zoom of the regulator dependence 
for dipole extrapolation based on CP-PACS data. Only the data 
point corresponding to the smallest value is shown at this scale. 



volume extrapolation such as shown in Figures |8] through 
IK)! The difference is that 100 infinite volume extrapola- 
tion points are produced close to the chiral regime. The 
exercise is to pretend that it is actual lattice QCD data. 
Clearly, a regularization scheme must be chosen to pro- 
duce the pseudodata. In this dipole regulator was 
chosen and pseudodata were created at A c = 1.0 GeV. 

The regularization dependence of the extrapolation is 
characterized by the scale dependence of the renormal- 
ized constants Cj. Consider how cq and Ci behave when 
analyzed with a variety of regulator values in Figures 1121 
and 1131 By choosing to use pseudodata produced at in- 
finite volume, one eliminates the concern that behaviour 
of the low energy constants across a range of regulators 
and pion masses is a finite volume artefact. The equiva- 
lents of Figures [T^] and Q2] for finite volume pseudodata 
exhibit the same features. 

If three pseudodata sets are compared, each with dif- 
ferent upper bounds on the range of m 2 . considered in 
the fit, an increasing regulator dependence in cq and ci 
is seen further from the PCR. A steep line indicates a 
strong scheme dependence in the result, and naturally 
occurs for data samples extending far outside the PCR. 
Scheme independence will appear as a completely hori- 
zontal graph. The latter is what one expects for a value 
of ml max < 0.04 GeV 2 , as shown in Figures IT2l and IT3l 
Note that, for each figure, all three curves (correspond- 
ing to different ml max ) arrive at stable values for cq and 
C2 on the righthand side of the graph (large A). To read 
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FIG. 12: (color online). Behaviour of cq vs. regulator parameter A, 
based on infinite volume pseudodata created with a dipole regulator 
at A c = 1.0 GeV (based on lightest four data points from PACS- 
CS). Each curve uses pseudodata with a different upper value of 
pion mass TO^ max . 
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FIG. 13: (color online). Behaviour of C2 vs. A, based on infinite 
volume pseudodata created with a dipole regulator at A c = 1.0 
GeV (based on lightest four data points from PACS-CS). Each 

curve uses pseudodata with a different upper value of pion mass 
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off the values of cq and C2 for large A is tempting but 
wrong. It is known that the correct values of cq and ci 
are recovered at A = 1.0 GeV, because at that value the 
pseudodata were created. 

The analysis of the pseudodata in Figures [T^] and [T3] 
shows that as the value of ml max is changed, the correct 
value of Co is recovered at exactly A = A c , where the 
curves intersect. This is also the intersection point for 
at A = A c . This suggests that when considering lattice 
QCD results extending outside the PCR, there may be 
an optimal finite-range cutoff. Physically, such a cutoff 
would be associated with an intrinsic scale reflecting the 
finite size of the source of the pion dressings. Mathemat- 
ically, this optimal cutoff is reflected by an independence 
of the fit parameters on ml max . 

To illustrate the non-triviality of this scale of curve- 
intersection, the pseudodata were analyzed with a differ- 
ent regulator, e.g. a triple-dipole regulator. Figures [T4l 
and [T51 show that the scale of the intersection is no longer 
a clear point, but a cluster centred about 0.5 to 0.6 GeV. 
The triple-dipole will of course predict a different 'best 
scale', since the shape of the regulator is different from 
that of the dipole used to create the pseudodata. The 
essential point of this exercise is that clustering of curve 
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FIG. 14: (color online). Behaviour of co vs. A, based on infinite 
volume pseudodata created with a dipole regulator at A c = 1.0 
GeV but subsequently analyzed using a triple-dipole regulator. 
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FIG. 15: (color online). Behaviour of C2 vs. A, based on infinite 
volume pseudodata created with a dipole regulator at A c = 1.0 
GeV but subsequently analyzed using a triple-dipole regulator. 



pion dressings of the hadron in question. 

The pseudodata analysis provides a good indication of 
a lower bound for A using a dipole regulator: Adi P > 
0.6 GeV. Similarly, Figure [T5l suggests a lower bound for 
the triple dipole regulator: A tr i p > 0.3 GeV. The same 
analysis can be repeated for the double dipole regulator 
to obtain A doub > 0.4 GeV. 

One may also constrain the lowest value that A can 
take by considering phenomenological arguments. Based 
on the physical values of the sigma commutator and the 
nucleon mass, a pion mass of ~ 0.5 GeV bounds 
the radius of convergence [H, H3, HU . Therefore, in or- 
der to ensure the inclusion of important contributions 
to the chiral physics, one should choose an energy scale 
A s harp ~ 0.5 GeV for a sharp cutoff (step function) reg- 
ulator. To compare this estimate for the sharp cutoff to 
that of dipole-like regulators, one can calculate the regu- 
lator value required such that u^(k 2 ) — 1/2 when the mo- 
mentum takes the energy scale of A s h ar p- This results in 
a rough estimate for a sensible value for the dipole, dou- 
ble dipole and triple dipole. These values are Adi P ~ 1.1 
GeV, Adoub ~ 0.76 GeV and A tr i p ~ 0.66 GeV, respec- 
tively. In any event, a wide range of regulator values will 
be considered, and the intersections of the curves for the 
low energy coefficients will be used in order to construct 
fits outside the PCR. This will be done in order to iden- 
tify the presence of an intrinsic scale for the pion source 
and an associated preferred regularization scale. 



intersections identifies a preferred renormalization scale 
that allows one to recover the correct low energy coeffi- 
cients. In this case, the crossing of the dash and dot-dash 
curves (from fitting) clearly identifies A^; p e = 0.6 GeV 
as a preferred regulator, which reflects the intrinsic scale 
used to create the data. Table U compares the values for 
Co and C2 recovered in this analysis for two different regu- 
lators: the preferred value A^ ip c = 0.6 GeV, and a large 
value Atri P = 2.4 GeV reflecting the asymptotic result 
recovered from DR. The input values of cq and ci used 
to create the pseudodata are also indicated. 

Note that the finite-range renormalization scheme 
breaks down if the finite-range regulator is too small. 
This is because A must be large enough to include the 
chiral physics being studied. The exact value of a sensi- 
ble lower bound in the finite-range regulator will depend 
on the functional form chosen as regulator. This is esti- 
mated for three dipole-like regulators in Section (IIV|) . 

Figure [T3l shows that the renormalization for C2 breaks 
down for small values of the regulator A. FRR breaks 
down for a value of Adi p much below 0.6 GeV, because the 
coefficients bt of the loop integral expansion in Eqs. (fl~8"|). 
(|2H)) and (|2"2"|) are proportional to A^ 3- ^. For high-order 
terms with large i, the coefficients will become large when 
A is small. This will adversely affect the convergence 
properties of the chiral expansion. One obtains a residual 
series expansion with good convergence properties only 
when A reflects the intrinsic scale of the source of the 



IV. INTRINSIC SCALE: LATTICE RESULTS 
A. Evidence for an Intrinsic Scale 

In the example of the pseudodata, an optimal finite- 
range cutoff was obtained from the data themselves. 
Clearly, the pseudodata have an intrinsic scale: the 
renormalization scale A c at which they were created. 
This test example leads the researcher to wonder if ac- 
tual lattice QCD data have an intrinsic cutoff scale em- 
bedded within them. That is, by analysing lattice QCD 
data in the same way as the pseudodata, can a similar 
intersection point be obtained from the renormalization 
scale flow of Co and C2? If so, this indicates that the lat- 
tice QCD data contain information regarding an optimal 
finite-range regularization scale, which can be calculated. 

The results for the renormalization of Co and ci as a 
function of A are now presented for JLQCD [25| , PACS- 
CS and CP-PACS [1| lattice QCD data. The 
JLQCD data use overlap fermions in two-flavor QCD, 
but the lattice box size for each data point is ~ 1.9 fm, 
smaller than the other two data sets. The PACS-CS data 
use the nonperturbatively 0(a)-improved Wilson quark 
action at a lattice box size of ~ 2.9 fm, but the data 
set only contains five data points and a large statistical 
error in the smallest m\ point. The CP-PACS data use 
a mean field improved clover quark action on lattice box 
sizes for each data point varying from ~ 2.2 fm to ~ 2.8 
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parameter 


input 


* Scale r\ r» 

Atrip =0.6 


Atrip =2.4 


Atrip = 2.4 






m ir,max = 0.25 


m 7T,max = 0.25 


m ir,max = 0.5 


CO 


0.902 


0.901 


0.899 


0.896 


C2 


3.00 


3.07 


3.17 


3.23 



TABLE I: A comparison of the parameters cq (GeV) and C2 (GeV ) at their input value (pseudodata created with a dipole at A c = 1.0 
GeV) with the values when analysed with a triple dipole regulator. Different values of A tr i p (GeV) and m 2 max (GeV 2 ) are chosen to 
demonstrate the scheme dependence of cq and C2 for data extending outside the PCR. 



fm. 

The chiral expansion is first used to chiral order 
0{m%). In this case, the fit parameters are Co and C2 
only. The results for a dipole regulator are shown in 
Figures [TBI through |2~T1 the results for the double dipole 
case are shown in Figures [2D through [57] and the results 
for the triple dipole are shown in Figures [28l through 1331 
To estimate the statistical error in the renormalized con- 
stants Sc, a bootstrap technique of 200 configurations of 
nucleon mass data is used. The configurations differ by 
the statistical error in the data, with values generated by 
a Gaussian distribution. In each plot, the same configu- 
rations are used for a variety of values of A considered. 
A few points are selected in Figures [16] through [33] to 
indicate the general size of the statistical error bars. 

It should be noted that none of these curves is flat 
to within 1% accuracy. All fits have included data be- 
yond the commonly accepted PCR. Clearly, there is a 
well-defined intersection point in the renormalization flow 
curves. Also, the value of A at which the intersection 
point occurs is the same even for different data sets, and 
for different c, . The tight groupings of the curve crossings 
lend credence to the ansatz of an intrinsic scale associ- 
ated with the finite size of the source of the pion dressings 
of the nucleon. This is a central result of this analysis. 

An intrinsic scale of A^ 10 w 1.3 GeV was obtained 
for the dipole, A| C Q ^ w 1.0 GeV for the double dipole 
and A^jp C w 0.9 GeV for the triple dipole. These values 
differ because the regulators have different shapes, and 
different values of A scalc are required to create a similar 
suppression of large loop momenta. 



B. Statistical Errors 

On each renormalization plot in Figures [TBI through |3"31 
there are many curves, each corresponding to different 
values of m^. max . It is of primary interest to what extent 
these curves match. Therefore, a Xdof should be con- 
structed, where dof equals the number of curves on each 
plot minus one for the best fit value of Co or C2, denoted 
by c T in the following. This also serves to quantify the 
constraint on the intrinsic scale A scale . The Xdof is eval- 
uated separately for each renormalized constant c (with 



error 5c) and regulator value A: 



1 - 

i=l 



(MA)) 2 



(31) 



for i corresponding to data sets with differing max . 
The theoretical value c T is given by the weighted mean: 



'(A) 



Eti^(A)/(MA)) 2 

ELi/(fc,(A)) 2 



(32) 



The Xdof can he calculated as a function of the regu- 
lator parameter A for each of the renormalization plots 
of Figures [TB] through OH In the case of the PACS-CS 
data, the minimum of the Xdof curve will be centred at 
the intersection point. In the case of the JLQCD and 
CP-PACS data, there appears to be a single intersection 
point on each plot, but in fact there are multiple intersec- 
tions over a very small window of A. The results for x\ Q f 
will indicate the 'best' central value of A. This central 
value of A will be taken to be the intrinsic scale. The 
Xdof curv es for a dipole regulator are shown in Figures 
[511 through |3"91 the x\ f curv es for the double dipole case 
are shown in Figures 00] through [55] and the x\ f curves 
for the triple dipole are shown in Figures 0BJ through 1511 



C. Higher Chiral Order 

Consider the renormalization of Co and C2 as a func- 
tion of A, for chiral order 0{m\ logm^). The results for 
PACS-CS and CP-PACS data are shown in Figures [5J 
through 1551 as an example. In this case, no clear inter- 
section points in the renormalization curves can be found, 
and so one is unable to specify an intrinsic scale. This cer- 
tainly should be the case when working with data entirely 
within the PCR, because all renormalization procedures 
would be equivalent (to a prescribed level of accuracy) 
and so there would be no optimal regulator parameter. 
It is known that this is not the case for the data sets used 
in this study. This is verified by considering the evident 
scale dependence of Co and C2 in Figures [55] through [55] 
The fact that Co and C2 change over the range of A values 
indicates that the data are not inside the PCR. Further, 
since no preferred scale is revealed, any choice of A ap- 
pears equivalent at this order. While this is encouraging 
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FIG. 16: (color online). Behaviour of cq vs. A, based on JLQCD 
data. The chiral expansion is taken to order O(m^), and a dipole 
regulator is used. A few points are selected to indicate the general 
size of the statistical error bars. 
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FIG. 20: (color online). Behaviour of cq vs. A, based on CP-PACS 
data. The chiral expansion is taken to order O(m^), and a dipole 
regulator is used. A few points are selected to indicate the general 
size of the statistical error bars. 
6.5 



5.5 - 




0.5 1 1 1 1 1 1 1 

0.0 0.4 0.8 1.2 1.6 2.0 2.4 
A (GeV) 

FIG. 17: (color online). Behaviour of C2 vs. A, based on JLQCD 
data. The chiral expansion is taken to order O(m^), and a dipole 
regulator is used. A few points are selected to indicate the general 
size of the statistical error bars. 
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FIG. 18: (color online). Behaviour of cq vs. A, based on PACS-CS 
data. The chiral expansion is taken to order 0(m^.), and a dipole 
regulator is used. A few points are selected to indicate the general 
size of the statistical error bars. 
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FIG. 19: (color online). Behaviour of C2 vs. A, based on PACS-CS 
data. The chiral expansion is taken to order O(m^), and a dipole 
regulator is used. A few points are selected to indicate the general 
size of the statistical error bars. 
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FIG. 21: (color online). Behaviour of C2 vs. A, based on CP-PACS 
data. The chiral expansion is taken to order O(m^), and a dipole 
regulator is used. A few points are selected to indicate the general 
size of the statistical error bars. 



that the scheme dependence is being weakened by work- 
ing to higher order, it must be recognized that there is a 
systematic error associated with the choice of A. In the 
case of the CP-PACS results shown in Figures [54] and 
[551 it can be seen that the statistical errors are substan- 
tially smaller than the systematic error associated with 
a characteristic range, Ai OW cr < A < oo, where Ai OW er is 
the lowest reasonable value of A. 

Since it is difficult to identify the intrinsic scale at this 
chiral order, the results for chiral order 0(m^) will be 
chosen to demonstrate the process of handling the exis- 
tence of an optimal regulator scale in lattice QCD data. 
The results for the calculation of the intrinsic scales A scalc 
for different data sets and regulators are given in Table 
HTl This table simply summarizes the central values from 
Figures [34] through 1511 Such excellent agreement be- 
tween the Co analysis and the C2 analysis is remarkable, 
and indicative of the existence of an instrinsic scale in 
the data. There is also consistency among independent 
data sets. It is important to realize that the value of 
A scalc is always the order of ~ 1 GeV, not 10 GeV, nor 
100 GeV; nor is it infinity. 

In calculating the systematic uncertainty in the ob- 
servables Co, C2 and the nucleon mass at the physical 
point due to the intrinsic scale at order O(m^logm ff ), 
two methods are provided. Firstly, the upper and lower 
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FIG. 22: (color online). Behaviour of cq vs. A, based on JLQCD 
data. The chiral expansion is taken to order 0(m^.) and a double 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 23: (color online). Behaviour of C2 vs. A, based on JLQCD 
data. The chiral expansion is taken to order 0(m^) and a double 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 24: (color online). Behaviour of cq vs. A, based on PACS-CS 
data. The chiral expansion is taken to order 0(m^.) and a double 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 26: (color online). Behaviour of cq vs. A, based on CP-PACS 
data. The chiral expansion is taken to order O(m^) and a double 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 27: (color online). Behaviour of C2 vs. A, based on CP-PACS 
data. The chiral expansion is taken to order O(m^) and a double 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 



regulator form 
optimal scale dipole double triple 



A scale 
cq , JLQCD 


1.44 


1.08 


0.96 


A scale 
C2 , JLQCD 


1.40 


1.05 


0.94 


A scale 
cq ,PACS — CS 


1.21 


0.93 


0.83 


A scale 

A c 2 , PACS-CS 


1.21 


0.93 


0.83 


A scale 

A c , CP-PACS 


1.20 


0.98 


0.88 


A scale 

A c 2 , CP-PACS 


1.19 


0.97 


0.87 



TABLE II: Values of the central A value in GeV, taken from the 
Xdof analysis for cq and C2, based on JLQCD, PACS-CS and CP- 
PACS data. 
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FIG. 25: (color online). Behaviour of C2 vs. A, based on PACS-CS 
data. The chiral expansion is taken to order 0(m^) and a double 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 



bounds from the Xdof analysis at order 0(m^) will be 
used to constrain A, and taken to be an accurate esti- 
mate of the systematic uncertainty in the contributions 
of higher order terms. Secondly, variation of the ob- 
servables across the characteristic range of scale values, 
Aiowcr < A < oo will be used, where Ai owor takes the value 
of 0.6, 0.4 and 0.3 GeV for the dipole, double dipole and 
triple dipole regulator, respectively. The results from 
both of these methods are displayed in Table IIIII 

The final results for the calculation of the renormalized 
constants cq, ci and the nucleon mass extrapolated to 
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FIG. 28: (color online). Behaviour of cq vs. A, based on JLQCD 
data. The chiral expansion is taken to order 0(mij.) and a triple 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 29: (color online). Behaviour of C2 vs. A, based on JLQCD 
data. The chiral expansion is taken to order 0(m^) and a triple 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 30: (color online). Behaviour of cq vs. A, based on PACS-CS 
data. The chiral expansion is taken to order 0(m^) and a triple 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 31: (color online). Behaviour of C2 vs. A, based on PACS-CS 
data. The chiral expansion is taken to order 0(m^) and a triple 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 32: (color online). Behaviour of cq vs. A, based on CP-PACS 
data. The chiral expansion is taken to order 0{m^) and a triple 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 33: (color online). Behaviour of C2 vs. A, based on CP-PACS 
data. The chiral expansion is taken to order 0(m^) and a triple 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 



the physical point (m^phys = 140 MeV) are summarized 
in Table IIV1 The lightest four data points from each 
of JLQCD, PACS-CS and CP-PACS lattice QCD data 
are used. The nucleon mass is calculated at the scale 
determined by the data. 



V. CONCLUSION 

In conclusion, it has been demonstrated that chiral ef- 
fective field theory is an important tool for investigat- 
ing the chiral properties of hadrons, and for extrapolat- 
ing lattice QCD results. Because the chiral expansion is 
only convergent within a PCR, a renormalization scheme 
such as finite-range regularization should be used for cur- 
rent lattice QCD results, and into the foreseeable future. 
Renormalization scheme dependence occurs when lattice 
QCD data extending outside the PCR are used in the 
extrapolation. This provides a new quantitative test for 
determining when lattice QCD data lie within the PCR. 
As most lattice data extend beyond the PCR, a formal- 
ism was developed to determine if there is an optimal 
regularization scale A scalc in the finite-range regulator, 
and to calculate it if it exists. It was concluded that such 
an optimal scale can be obtained from the data itself by 
analyzing the renormalization flow curves of the low en- 
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regulator form 




sys. err. 


dipole 


double 


triple 



rA JLQCD 
c 





.001, 





.009 





.001, 





.013 





.001, 





.016 


cA„PACS-CS 

o c 





.005, 





.006 





.005, 





.010 





.006, 





.012 


rA CP — PACS 

o c 





.002, 





.024 





.002, 





.037 





.002, 





.045 


rA JLQCD 
c 2 





.02, 





.31 





.03, 





.38 





.01, 





.48 


A A„PACS-CS 
o c 2 





.18, 





.25 





.16, 





.33 





.14, 





.43 


rA„CP-PACS 





.02, 





.40 





.02, 





.58 





.02, 





.73 


lA M JLQCD 
° JU iV,phys 





.0004, 





.0051 





.0003, 





.0073 





.0003, 





.0090 


rA.^PACS-CS 
JU iV,phy S 





.0022, 





.0030 





.0025, 





.0046 





.0025, 





.0058 


tAifCP-PACS 
5 M IV,phys 





.0012, 





.0175 





.0013, 





.0270 





.0014, 





.0326 



TABLE III: Results at 0(ra\ logm,^) for the systematic error due to the intrinsic scale, calculated using two methods, for the values of 
co (GeV), C2 (GeV -1 ) and the nucleon mass Mjv (GeV) extrapolated to the physical point (m^ , p hys = 140 MeV). The first number in 
each column is the systematic error due to the intrinsic scale using the upper and lower bound from the x\ a j analysis at order 0(m^). 
The second number is the systematic error due to the instrinsic scale across the whole range of A values from the lowest reasonable value 
(A = Ai ower ) obtained from the pseudodata analysis, to the asymptotic value (A = co). 







regulator form 






parameter 


dipole 


double 


triple 


WM(1) 


WM(2) 


JLQCD 
C 


0.873(18)(16) 


0.875(17)(16) 


0.891(17)(16) 


0.880(29) 


0.879(32) 


„PACS-CS 
C 


0.900(51)(15) 


0.899(51)(14) 


0.898(51)(14) 


0.899(53) 


0.899(55) 


„CP-PACS 
C 


0.924(3)(8) 


0.914(3)(7) 


0.918(3)(7) 


0.918(13) 


0.920(37) 


JLQCD 
C 2 


3.09(9)(11) 


3.18(9)(12) 


3.20(9)(11) 


3.16(18) 


3.14(43) 


„PACS-CS 
c 2 


3.06(32)(15) 


3.15(31)(14) 


3.17(31)(14) 


3.13(39) 


3.12(49) 


„CP-PACS 
c 2 


2.54(5) (4) 


2.70(5)(2) 


2.71(5)(3) 


2.66(18) 


2.61(60) 


a j-JLQCD 

JW JV,phys 


1.02(2)(9) 


1.02(2)(9) 


1.02(2)(9) 


1.02(9) 


1.02(9) 


a i-PACS-CS 
IV1 N,phys 


0.967(45) (43) 


0.966(45) (43) 


0.966(45) (43) 


0.966(62) 


0.966(62) 


» j-CP-PACS 
1V1 N,phys 


0.982(2) (40) 


0.975(2)(43) 


0.978(2)(42) 


0.979(43) 


0.979(50) 



TABLE IV: Results at 0(m% log m n ) for the values of co (GeV), C2 (GeV 1 ) and the nucleon mass Mjv (GeV) extrapolated to the 
physical point (rrt^phys = 140 MeV). WM is the weighted mean of each row. The nucleon mass is calculated at the optimal scale A scale , 
which is the average of A^ ale and A^ ale for each data set. The extrapolations are performed at box sizes relevant to each data set: 

^cxtrap 3 = L9 fm > L rrf?a S p ~ CS = 2 - 9 fm and L £arlip ACS = 2 - 8 fm - The errors are quoted as the estimate of the statistical error first 
(based on random bootstrap configurations), and the systematic error obtained from the number of m^. values used second. Two seperate 
weighted means are calculated for each row. WM(1) incorporates the systematic error in the intrinsic scale using the upper and lower 
bound from the Xdof ana b/ s i s a t order O(m^). The WM(2) incorporates the systematic error due to the intrinsic scale across the whole 
range of A values from the lowest reasonable value (A = Ai owcr ) obtained from the pseudodata analysis, to the asymptotic value (A = oo). 
The weighted means also include an estimate of the systematic error in the choice of regulator. All errors are added in quadrature. Note 
that any order 0(a) errors have not been incorporated into the total error analysis. 



ergy coefficients in the chiral expansion. The optimal 
scale is selected by the value for which the renormalized 
constants are independent of the upper bound of the fit 
domain. This also means that the renormalized constants 
are not to be identified with their asymptotic values at 
large A. 

It was revealed that a preferred regularization scheme 
exists only for data sets extending outside the PCR. Such 
a preferred regularization scheme is associated with an 
intrinsic scale for the size of the pion dressings of the 
nucleon. By working to sufficiently high chiral order, it 
was discovered that the scale dependence was weakened. 



Nevertheless, the residual scale dependence persists as a 
significant component of the systematic uncertainty. For 
efficient propagation of this uncertainty, an interesting fu- 
ture direction would be to consider marginalization over 
the scale dependence [2^|. The described procedure was 
used to calculate the nucleon mass at the physical point, 
the low energy coefficients cq and C2 and their associ- 
ated statistical and systematic errors. Several different 
functional forms of regulator were considered, and lattice 
QCD data from JLQCD, PACS-CS and CP-PACS were 
used. An optimal cutoff scale A scale for each set of lattice 
QCD data was obtained, and the systematic error in the 
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FIG. 34: (color online). Behaviour of Xdof ^ or c ° VSi ^> based on 
JLQCD data. The chiral expansion is taken to order 0(m^.), and 
a dipolc regulator is used. 
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FIG. 38: (color online). Behaviour of Xdof ^ or c ° VSi ^> Dase d 
on CP-PACS data. The chiral expansion is taken to order 0(m'f r ), 
and a dipole regulator is used. 
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FIG. 35: (color online). Behaviour of Xdof ^ or C2 vs ' ^> Dase d on 
JLQCD data. The chiral expansion is taken to order O(m^), and 
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FIG. 39: (color online). Behaviour of x\ a f ^ or C2 vs ' ^> Dase d 
on CP-PACS data. The chiral expansion is taken to order 0(m^), 
and a dipole regulator is used. 
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FIG. 36: (color online). Behaviour of x\ a f f° r °0 vs. A, based on 
PACS-CS data. The chiral expansion is taken to order O(m^), and 
a dipole regulator is used. 
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FIG. 37: (color online). Behaviour of x\ f ^ or C2 VSi ^» based on 
PACS-CS data. The chiral expansion is taken to order O(m^), and 
a dipole regulator is used. 



In summary, the existence of a well defined intrinsic 
scale has been discovered. It has also been illustrated how 
its value can be determined from lattice QCD results. 
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FIG. 40: (color online). Behaviour of Xdof ^ or c ° VSi ^> based on 
JLQCD data. The chiral expansion is taken to order 0(m^.), and 
a double dipole regulator is used. 
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FIG. 41: (color online). Behaviour of Xdof f° r c 2 vs. A, based on 
JLQCD data. The chiral expansion is taken to order O(m^), and 
a double dipole regulator is used. 
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FIG. 42: (color online). Behaviour of Xdof ^ or c ° vs - A, based on 
PACS-CS data. The chiral expansion is taken to order O(m^), and 
a double dipole regulator is used. 
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FIG. 43: (color online). Behaviour of Xdof ^ or C2 VSi ^» based on 
PACS-CS data. The chiral expansion is taken to order O(m^), and 
a double dipole regulator is used. 
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FIG. 44: (color online). Behaviour of Xdof ^ or c ° VSi ^> based 
on CP-PACS data. The chiral expansion is taken to order 0(m^.), 
and a double dipole regulator is used. 
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FIG. 45: (color online). Behaviour of x\ j ^ or ° 2 vs ' ^> based 
on CP-PACS data. The chiral expansion is taken to order 0(m^.), 
and a double dipole regulator is used. 
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FIG. 46: (color online). Behaviour of x\ a f ^ or c ° VSi ^> Dase cl on 
JLQCD data. The chiral expansion is taken to order 0(m^.), and 
a triple dipole regulator is used. 
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FIG. 50: (color online). Behaviour of Xdof ^ or c ° VSi ^> Dase d 
on CP-PACS data. The chiral expansion is taken to order 0(m'f r ), 
and a triple dipole regulator is used. 
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FIG. 47: (color online). Behaviour of Xdof ^ or C2 vs ' ^> Dase( l on 
JLQCD data. The chiral expansion is taken to order O(m^), and 
a triple dipole regulator is used. 
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FIG. 51: (color online). Behaviour of x\ a f ^ or C2 vs ' ^> Dase< l 
on CP-PACS data. The chiral expansion is taken to order 0(m^), 
and a triple dipole regulator is used. 
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(color online). Behaviour of Xdof ^ or c ° vs ' ^> Dase£ i on 
PACS-CS data. The chiral expansion is taken to order O(m^), and 
a triple dipole regulator is used. 
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FIG. 49: (color online). Behaviour of x\ f ^ or C2 VSi ^» bas e d on 
PACS-CS data. The chiral expansion is taken to order O(m^), and 
a triple dipole regulator is used. 
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FIG. 52: (color online). Behaviour of cq vs. A, based on PACS-CS 
data. The chiral expansion taken to order 0(m*logm7r) and a 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 53: (color online). Behaviour of C2 vs. A, based on PACS-CS 
data. The chiral expansion taken to order 0(m^ log m^) and a 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 54: (color online). Behaviour of cq vs. A, based on CP-PACS 
data. The chiral expansion taken to order 0(m* logm,,-) and a 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 
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FIG. 55: (color online), 
data. The chiral expansion taken to order 0(m^ log m^) and a 
dipole regulator is used. A few points are selected to indicate the 
general size of the statistical error bars. 



